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Abstract. The cosmic microwave background radiation is supposed to be Gaussian and this 
hypothesis is in good agreement with the recent very accurate measurements. Nonetheless a 
tiny amount of non-Gaussianity is predicted by the standard inflation scenario, while more 
exotic models suggest a higher degree of non-Gaussianity. Tightly constraining the level 
of Gaussianity in the CMB data represents then a fundamental handle to understand the 
physics and the origin of our universe. By means of needlets, a novel rendition of wavelets, 
characterised by excellent properties of localisations both in harmonic and pixel domain, we 
are able to detect anomalous spots in the southern hemisphere responsible for roughly the 
50% of power asymmetry we measure in the CMB power spectrum, and to perform a de- 
tailed analysis of the needlets bispectrum. We then constrain the primordial non-Gaussianity 
parameter, /nl = 21 ± 40 at 68% c.f., and spot a high asymmetry in the bispectrum, in par- 
ticular in the isosceles configurations. 
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1. Introduction 

During the last decade the quality of the cos- 
mological observations has increased impres- 
sively, providing theorists with very accurate 
datasets. Cosmic Microwave Background ra- 
diat ion (CMB) measur ements j Hinshaw et al.l 
|2009), SuperNovae lA (Riess et al. 2009) and 
Baryo nic Acoustic Oscillations ( Percival et al.i 
|2009|) fit pretty well within the scenario of the 
so-called Cosmological Concordance Model 
(ACDM), which describes fairly well the evo- 
lution of the Universe by means of an handful 
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of parameters (iKomatsu et al.ll2009l) . We live 
in a flat universe, whose critical density is pro- 
vided by baryons Qb - 4%, dark matter Q.^ - 
25%, responsible for the structures formation, 
and the troublesome vacuum energy, Qa - 
70%. The other parameters set the normalisa- 
tion. As, and the power, n^, of the primordial 
cosmological fluctuations and the contribution 
of the stars formation at more recent epoch, t. 
Despite its simplicity, the ACDM model lacks 
of a solid theoretical basis. The value and the 
origin of the v acuum energy is far from be- 
ing understood jCopeland et al.ll2006h . and the 
mechanism which seeds the cosmological per- 
turbations, namely inflation (iGutli,1981i) . re- 
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mains more a scenario than a full tested the- 
ory. Our Universe is supposed to be homo- 
geneous and isotropic with nearly Gaussian 
fluctuations. Indeed this is in good agreement 
with the most recent data, even though the 
mechanism generating the perturbations itself 
should intro duce a small no n-Gaussian con- 
tribution (Ba rtolo et a n 12004 . Moreover, ex- 
otic early univers e scenarios, such as brane- 
inspired m odels ( Stein hardt & TuroS |2002[) . 
multi-fie lds (iLvth & Wa nds 2002) or ekpyrotic 
models jMizuno et alj 2008) predict a higher 
level of non-Gaussianity. Accurately determin- 
ing the statistics of the cosmological pertur- 
bations represents a fundamental handle on 
the early universe physics necessary for under- 
standing the nature and the evolution of our 
Universe. 

In the following we tackle the issue of 
non-Gaussianity in the CMB data, analysing 
the WMAP 5 -year temperature data by means 
of needlets. Needlets are a novel rendition of 
wa velets introduced first in functional analy- 
sis (iNarcowich et al. 2006) and then extended 
to statistical (Baldi et al. 2006) and CMB data 
analysis ( Pietrobon et al.. 2006 ). Needlets are a 
basis (more properly a tight frame) defined di- 
rectly on the sphere, which shows an exponen- 
tial localisation property both in pixel and har- 
monic space. Moreover they are very weakly 
correlated: this makes them particularly suit- 
able for CMB data analysis where partial sky 
coverage, noise and beam effects have to be 
minimised in order to extract tiny signals, es- 
pecially when non-Gaussian. Needlets are a 
quadratic combination of spherical harmonics, 
weighted by a window function, b{ilB^) 



with the needlet profile in real space. The local- 
isation property is clearly visible. 
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Fig. 1. Profile of the b{e/Bj) filter (top) and needlet 
(bottom) for the choice B = 1.8 j = 4. 



2. Spot detection and Power 
spectrum asymmetry 

Since the first release, the WMAP data 
have been tested ag ainst Gaussian i ty ari d 
asymme try (see e.g. 'Eriksen et al.' ('2004*); 
Cruz et a l. (2005); Land & Magueijo (2007); 
de Oli veira-Costa et al.l (I2004i) ). We applied 
needlets to the 5-year data to study the 
anomaly/asymmetry problem in a coherent 
framework. We extracted the needlets coeffi- 
cients from the temperature map given by 



where y represents a direction in the sky, {^jk} 
is a set of cubature points (for practical pur- 
pose iden tified with a pixelization of the sphere 
dGorski et al.ll2005h l and B a user chosen pa- 
rameter which determines the filter width in 
harmonic space . For a detailed discussion see 
iMarinucci et alj (I2008h and refs. therein. An 
example of this filter is given in Fig.[T] together 



Ik)- 



(2) 



Pjk can be easily represented in a moU- 
weide projection: the case of the Internal 
Linear Combination maf0 is given in Fig. |2l 
Notice that needlets coefficients are analyti- 



http://lambda.gsfc.nasa.gOv/~/ilc_map_get.cfm 



Pietrobon: NG seen through Needlets 



3 



Fig. 2. Needlet coefficients of the WMAP 5-year 
CMB temperature map. 6 = 1.8 and j = 4. 



cally linked to those of the spherical harmon- 
ics: this guarantees an expression for the 2- 
and 3-point correlation functions, respectively 
the power spectrum and bispectrum. By ap- 
plying needlets, we found three very signif- 
icant spots in the southern hemisphere (the 
brightest features in Fig. |2l) which seem to be 
bar ely compatible with th e Gaussian hypothe- 
sis dPietrobon et al.ll2008h . Thanks to needlets 
localisation, we were able to determine which 
angular scales they span. In particular comput- 
ing the needlets power spectrum 

where Q - (|afmP), we found that the sig- 
nificant power asymmetry between the north 
and the south hemisphere is localised at large 
scales, where the spots peak. Moreover, mask- 
ing the spots, we measured a decrease of the 
difference in power of a factor 2, underlin- 
ing the substantial contribution given by these 
anomalous features. Fig. |3] summarises our re- 
sults. We checked whether this effect influ- 
ences the parameter estimation, finding that, 
with the current precision, the difference re- 
mains within one sigma. 



3. Bispectrum analysis 

The signature of non-Gaussianity appears in 
the higher moments of a distribution, which 
are no longer completely specified by the first 
moment (i.e. the mean value of the distribu- 
tion) and the second moment (i.e. the stan- 
dard deviation). For a Gaussian distribution. 
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Fig. 2).Pj extracted from WMAP data for different 
sky cuts. The black solid line shows the signal for 
the whole CMB sky (kq85 is applied). The blue and 
red dashed lines show how the power is split be- 
tween the two hemispheres. When the cold and hot 
spots detected are masked the excess of power in the 
southern region is decreased (light blue and orange 
dot-dashed lines). 

all odd moments are vanishing, while the even 
ones can be expressed in term of just the first 
two. We then look for a non-vanishing bispec- 
trum of the distribution of the needlets coeffi- 
cients. The amplitude of the bispectrum is usu- 
ally parameterised by the non-linear parameter 
/nl which governs the amplitude of the non- 
Gaussian contribution to the primordial gravi- 
tational potential expansion with respect to the 
linear leading bit: 

<t(x) = Ol(^) + /nl [DlW - ■ (4) 

The primordial fluctuations are converted into 
CMB ones, whose bispectrum is given by 
{a{^,n^at,j„^_a(^m^). In needlets space the former 
expression reads: 



(5) 



where A^pix is the number of pixels outside 
the applied mask and crj the variance of the 
needlets coeffic ient at the j resolution. See 
iPietrobon et al.l ([2009.) and Refs. therein for 
a detailed discussion. We extracted needlets 
coefficients from the noise weighted combi- 
nation of the WMAP5 channels and mea- 
sured fNi, by appl ying the following estimator 
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/nl - 



(6) 
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We obtained /nl = 21+40 at 68% c.l. We cal- 
ibrated our result by means of b oth Gaussian 
and non-Gaussian ( iLiguori et al.l2 007) simula- 
tions of the underlying ACDM best fit model. 
Our constraints are consistent with those ob- 
tained with different techniques: having several 
tools, affected by different systematics, is cru- 
cial for analysing the upcoming cosmological 
experiments datasets. 

The /nl parameter encodes all the informa- 
tion contained in the bispectrum irrespective of 
the specific triangle configuration. The bispec- 
trum signal can be actually split according to 
the geometry in equilateral, isosceles, scalene 
and open configurations. Different kinds of 
non-Gaussianity may result in different shapes, 
so it is worth looking at them s eparately. 
This i s what we addressed in Pietrobon et alJ 
(120091) . The main result of the paper is sum- 
marised in Tab. [3] where the percentage of ran- 
dom Gaussian simulations with a higher 
than the data is quoted. Again, the southern 



conf. FULL SKY NORTH SOUTH 

all (115) 29% 96% 2% 

equi(9) 20% 11% 45% 

iso (56) 5% 96% 0.5% 

seal (50) 60% 90% 7% 

open (50) 3% 85% 2% 



Table 1. Percentage of the simulations with a 
larger than WMAP 5-year data for the different tri- 
angular configurations of the needlets bispectrum. 



hemisphere results very anomalous, especially 
in the isosceles configurations, those charac- 
terised by a local type of non-Gaussianity. 

4. Conclusions 

We have discussed the importance of char- 
acterising the statistical distribution of the 
cosmological perturbations to understand the 
early universe physics. We study in detail 
several properties of the CMB sky where 
a non-Gaussian signal may arise by means 
of needlets. We found anomalous spots in 
the southern hemisphere which account for 



the 50% of the power spectrum asymme- 
tries. Moreover we measured the needlets 
bispectrum constraining the primordial non- 
Gaussianity parameter and studying its prop- 
erties according to the geometry of the triangle 
configurations. The isosceles ones turn out to 
be the most anomalous: whether this is a sig- 
nature of a peculiar early universe model is an 
interesting issue in light of the new upcoming 
experiments. 
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